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We demonstrate theoretical conditions for highly-efficient degenerate four-wave mixing in triply- 
resonant nonlinear (Kerr) cavities. We employ a general and accurate temporal coupled-mode 
analysis in which the interaction of light in arbitrary microcavities is expressed in terms a set of 
coupling coefficients that we rigorously derive from the full Maxwell equations. Using the coupled- 
mode theory, we show that light consisting of an input signal of frequency wq — Aw can, in the 
presence of pump light at wo, be converted with quantum-limited efficiency into an output shifted 
signal of frequency wo + Aw, and we derive expressions for the critical input powers at which 
this occurs. We find that critical powers in the order of lOmW assuming very conservative cavity 
parameters (modal volumes ~ 10 cubic wavelengths and quality factors ~ 1000. The standard 
Manley-Rowe efficiency limits are obtained from the solution of the classical coupled-mode equations, 
although we also derive them from simple photon-counting "quantum" arguments. Finally, using a 
linear stability analysis, we demonstrate that maximal conversion efficiency can be retained even in 
the presence of self- and cross-phase modulation effects that generally act to disrupt the resonance 
condition. 

PACS numbers: 42.65.Ky, 42.60.Da, 42.65.Sf, 42.65.Jx 



I. INTRODUCTION 

For many years, researchers have used confinement of 
light for a long time in a small volume (resonant cavities) 
to enhance light-matter interactions such as optical non- 
linearities, recently entering the integrated-optics regime 
of smaller and smaller cavities with limited sets of inter- 
acting modes. In such systems, careful design is required 
to maximize the efficiency and minimize the power of a a 
given nonlinear process such as frequency conversion [l|- 
[l7| , and the use of cavities can also lead to qualitati vely 
new phenomena such as bi/multistability [l], 0, [l8l43l| . 
While the use of cavities is known to enhance nonlinear 
effects, every distinct nonlinear process requires a new 
analysis. In this paper, we consider the problem of intra- 
cavity degenerate four- wave mixing (DFWM): an elec- 
tromagnetic cavity resonant at three frequencies wo and 
u}q ± Aw, in which a third-order (x ) nonlinearity con- 
verts an input signal at luq — Aw to an output shifted 
signal at loq + Aw in the presence of an input pump 
at ujq. The small- Aw regime corresponds, for example, 
to conversion between different channels in wavelength- 
division multiplexing (WDM), similar to recent experi- 
mental studies of nonlinear frequency conversion in cav- 
ities [32h4{)| ; it allows one to exploit structures like ring 
resonators[3| or photonic-crystal cavities [4l| that sup- 
port closely spaced resonances, and in Sec. II V Al we show 
that this regime supports stable quantum-limited con- 
version at low signal powers for a critical pump power. 
Conversely, we show that the Aw > wo regime gener- 
alizes our previous work on intra-cavity third-harmonic 
generation (THG) [J 0, and in Sec. llVBl we show that 
this regime supports stable conversion with 100% effi- 
ciency at a critical pump and signal power. For exam- 



ple, with a typical nonlinear material such as gallium 
arsenide (GaAs) and reasonable cavity parameters (vol- 
ume ~ 10 cubic wavelengths and quality factors ~ 1000), 
we obtain critical powers in the milliwatts (on the order 
of lOmW) for both Aw regimes. The standard Manley- 
Rowe efficiency limits are considered from both a sim- 
ple photon-counting "quantum" argument flU l43j and 
are also derived from purely classical coupled-mode equa- 
tions 0, EE} (Sec. EI)) , where the latter also yield sta- 
bility information, critical powers, and other dynamics 
(Sees. ITV^l and HV Bp . The coupling coefficients in these 
equations are derived explicitly from the full Maxwell 
equations for arbitrary microcavities (Sec. III]) . We also 
show that the nonlinear dynamics lead to additional 
phenomena, such as multistability and limit-cycle (self- 
pulsing) solutions, similar to phenomena that were pre- 
viously shown for other nonlinear systems 0, EH, S3 
(Sees. ITVXl and ITVFl)) . Finally, in Sec. lTVCl we consider 
the effects of self- and cross-phase modulation (SPM and 
XPM), which induce nonlinear shifts in the cavity fre- 
quencies: these must be compensated by pre-shifting the 
resonances and also affect the stability analysis (as we 
previously found for THG [l[). 

Although nonlinear effects in electromagnetism are 
weak, it is well known that confining light in a small vol- 
ume and/or for a long time, as in a waveguide or cavity, 
can both enhance the strength and modify the nature 
of nonlinear phenomena (48l |49| . Much previous work 
in nonlinear frequency conversion has studied x P r0 ~ 
cesses (where there is a change in the susceptibility that 
is proportional to the square of the electric field) such as 
second harmonic generation (SHG) I21 13. IItI lia l5(il - [59j . 
difference-frequency generation (DFG) Oa, la l60W62l ] . and 
optical parametric amplification (OPA) [63146 5] . Stud- 
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FIG. 1: (Color online) (Left) Schematic for degenerate four- 
wave mixing involving a coupled waveguide-cavity system. 
Dynamical variables for coupled-mode equations represent: a 
single input (output) channel (with incoming/outgoing field 
amplitudes s±) coupled to a resonant cavity with three modes 
at frequencies uio, u) m = cjo — Auj and ui p — ujq + Aui (and 
corresponding amplitudes do, a m and a p ). The three reso- 
nant modes are nonlinearly coupled by a Kerr (x' 3 ') nonlin- 
earity. (Right) Diagram illustrating the relationship between 
the three resonant frequencies. 



ics of SHG in doubly resonant x*' 2 ' cavities have demon- 
strated that 100% conversion efficiency is achieved at 
critical pump power, much lower than for SHG in singly 
resonant cavities @, [3, H3, HI, H^-E!] . Recent stud- 

ies of DFG in triply resonant \^ cavities also showed 
the existence of a critical relationship between pump 
and idler power that results in optimal quantum-limited 
with potential applications to terahertz 
The existence of quantum-limited frc- 



conversion 
generation 

quency conversion can be predicted from the Manley- 
Rowe relations, which govern the rates of energy trans- 
fer in nonlinear systems |48j . There has also been some 
recent work on intra-cavity x^ third harmonic gener- 
ation (In a x*- 3 ) medium, there is a change in the 
refractive index proportional to the square of the elec- 
tric field.) As in SHG, THG in doubly resonant cavities 
has been shown to support solutions with 100% conver- 
sion efficiency, even when taking into account nonlinear 
frequency-shifting due to SPM and XPM, as well as in- 
teresting dynamical behavior such as multistability and 
limit cycles (self-pulsing) [l|, with lower power require- 
ments compared to singly resonant cavities or nonres- 
onant structures 0, [SOL l59l l65l . Ir39l l69j - [8l| . Limit cy- 
cles have been observed in a number of other nonlin- 
ear optical systems, including doubly resonant x*- 2 ' cavi- 
ties 0, , bistable multimode Kerr cavities with time- 
delayed nonlinearities [83fl. nonresonant distributed feed- 
back in Braggs gratings |l9[, and a number of nonlinear 
lasing devices [84j ■ 

In what follows, we extend the previous work on SHG, 
DFG and THG in resonant cavities to the case of DFWM 
in x'' 3 ' media. Four-wave mixing is characterized by tak- 
ing input light at frequencies Wi,<J2, and W3 and pro- 
ducing light at frequency W4 = ±wi ±^±1^3; degener- 
ate four-wave mixing, however, is restricted to the case 
where uj\ = UJ3 to generate 2wi — uj^. Previous work 
has studied FWM in the context of optical fibers [85l — 
[89[ and even matter waves [§(| , as well as demonstrating 
the use of FWM in applications such as phase conjuga- 
tion 6 ^ 91-93 | and generation of two-photon coherent 
states |85l. |94|. |95|] . While there has been recent exper- 



imental work on intra-cavity FWM in x^ media (de- 
generate or otherwise) [32T - |40l [96|-[98| , we are not aware 
of any detailed studies of the underlying theoretical phe- 
nomena in general cavities. As we shall, see DFWM in 
triply-resonant cavities shares many qualitative features 
with SHG, DFG, and THG, including the existence of 
critical powers at which optimal conversion efficiency is 
achieved as well as interesting nonlinear phenomena such 
as limit cycles and multistability. As in DFG, and unlike 
SHG or THG, there exist Manley-Rowe limitations on 
the overall conversion efficiency. In Sec. IIII1 we discuss 
the corresponding relations governing four-wave mixing 
and illustrate their implications for conversion efficiency. 
These relations can be obtained classically through tem- 
poral coupled-mode theory 0, 5H , but they are more 
easily motivated and understood from a quantum per- 
spective [13, HH . Such argumen ts ha ve been emp loyed 
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before in the context of lasing [99j, llOOj , RF circuits 
and other nonlinear optics phenomena [48|. In the case 
of intra-cavity frequency conversion, we show how both 
perspectives yield limits on conversion efficiency. 

Several different approaches can be used to study non- 
linear optical systems. Most directly, brute-force numer- 
ical simulations by a variety of met hods , such as finite- 
difference time-domain (FDTD) (ZllllO^. offer the most 
general and flexible technique, in that they can charac- 
terize phenomena involving many degrees of freedom and 
going beyond the pcrturbative regime, but such simu- 
lations are relatively slow and allow one to study only 
a single geometry and excitation at a time. More ab- 
stract analyses are possible in many problems because 
confinement to a waveguide or cavity limits the degrees 
of freedom to the amplitudes of a small set of normal 
modes, combined with the fact that optical nonlineari- 
ties are typically weak (so that they can be treated as 
small perturbations to the linear modes). For example, 
many nonlinear phenomena have been studied in the con- 
text of co-propagating plane waves, in which the am- 
plitudes of the waves can be shown to satisfy a set of 
simple ordinary differential equations (ODEs) in space 
(the slowly varying envelope approximation) [48j. More 
generally, however, it can be shown that all nonlinear 
problems coupling a finite set of modes and satisfying 
certain fundamental principles such as conservation of 
energy, regardless of the underlying wave equation (e.g. 
electromagnetic or acoustic waves), can be described by 
a universal set of ODEs characterized by a small number 
of coefficients, determined by the specific geometry and 
physics. This approach, which has come to be known 
as temporal co upled-mo de theory (TCMT), dates back 
several decades [44], Il03j | and has been applied to a large 
numb er o f problems, from microwave transmission sys- 
tems [l0l( | to the nonlinear intra-cavity prob lems (SHG, 
DFG, and THG) mentioned above [Ell M, GUI. Like- 
wise, we employ TCMT in this paper to characterize the 
most general possible behavior of intra-cavity DFWM 
systems, regardless of the nature of the cavity. As re- 
viewed elsewhere [44[ , TCMT begins with the purely lin- 
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ear system and breaks it into abstract components such 
as input/output channels (e.g. waveguides or external 
losses) and cavities, characterized by resonant frequen- 
cies and coupling rates that depend on the geometry; 
it then turns out that the ODEs describing such a sys- 
tem are completely determined by those parameters once 
the constraints of conservation of energy, linearity, time- 
invariance, and reciprocity (or time-reversal invariance) 
arc included, under the key assumption that coupling 
rates are slow compared to the frequencies (i.e., strong 
confinement) [U |4J] . Nonlinearities can then be intro- 
duced as additional terms in these equations, without 
disturbing the previously derived relationships, as long 
as the nonlinear processes are also weak (i.e. nonlinear 
effects occur slowly compared to the frequency), which 



is true in nonlinear optics |48[. Using these ODEs, the 



general possible behaviors can be obtained (including the 
Manlcy-Rowe relations mentioned above), but to obtain 
the specific characteristics of a particular geometry one 
then needs a separate calculation to obtain the cavity 
parameters. Properties of the linear modes such as fre- 
quencies and lifetimes (Q) can be obtained by standard 
computational methods [u], [44[. It turns out that the 
nonlinear coefficients can also be obtained from the lin- 
ear calculations, thanks to the fact that the nonlinearities 
are weak: using perturbation theory, expressions for the 
nonlinear coefficients as integrals of the linear modes can 
be derived from Maxwell's equations. Such expressions 
were previously derived for SHG and THG 0], and also 
recently for DFG [6(|. Here, we derive both the abstract 
TCMT equations and the specific nonlinear coupling co- 
efficients for DFWM in the Maxwell equations with 
nonlinearities. 

We begin apply the coupled-mode formalism to the 
case of DFWM in a triply-resonant cavity in Sec. [Til 
to obtain the coupled-mode equations of motion as well 
as explicit expressions for the nonlinear coupling coeffi- 
cients. We then briefly discuss general properties of the 
conversion process in Sec. IIIII and, using the standard 
Manlcy-Rowe relations and simple photon-counting ar- 
guments, obtain limits on the maximal efficiency of the 
system. In Sees. HVAl and [TV B[ we analyze the stabil- 
ity and dynamics of the solutions to the coupled-mode 
equations obtained in Sec. HH neglecting SPM and XPM 
effects, and demonstrate the existence of the maximal 
conversion efficiencies obtained in Sec. IIIII Finally, in 
Sec. IIVCI we briefly consider the effects of SPM and 
XPM using a simple model to illustrate the qualitative 
behavior of the system; in particular, wc demonstrate 
the existence of stable, maximal efficiency solutions even 
including SPM and XPM effects. 



II. TEMPORAL COUPLED-MODE THEORY 

We consider the situation depicted schematically in 
Fig. [1] an input/output channel coupled to a triply- 
resonant nonlinear (x' 3 ^) cavity. Here, input light at 



wo and ui m = loq — Aw is converted to output light at 
a different frequency w p = ujq + Aw, where Aw deter- 
mines the separation between the three frequencies. The 
frequency-conversion process occurs inside the nonlinear 
cavity, which supports resonant modes of frequencies wo, 
w m , and w p , and corresponding modal lifetimes (or 
quality factors Qk = WfcTfc/2 [44[) describing the overall 
decay rate (I/t^) of the modes. In particular, the to- 
tal decay rate consists of decay into the output channel, 
with rate l/T St k, as well as external losses (e.g. absorp- 
tion) with rate l/r e .fe, so that 1/tu = 1/ T *,fc + \/r e ^- 
Note that, to compensate for the effects of SPM/XPM, 
as described in [l[ and in Sec. II V CI we will eventually use 
slightly different cavity frequencies w£ av that have been 
prc-shiftcd away from W&. 

It is most convenient to express the TCMT equations 
in terms of the following degrees of freedom [u], |44| : we 
let etfc denote the time-dependent complex amplitude of 
the fcth mode, normalized so that |afc| 2 is the electromag- 
netic energy stored in this mode, and Sk,± denote the 
time-dependent amplitude of the incoming (+) and out- 
going (— ) wave, normalized so that |sfc.±| 2 is the power 
in the fcth mode. (In what follows, we take s p , + = 0, 
corresponding to the up-conversion of light at wo and w m 
to light at u> p) for Aw > 0. In order to study the alterna- 
tive down-conversion process, one has but to set Aw < 0, 
in which case we effectively have w p — > w m , as described 
below.) 

The derivation of the linear TCMT equations, corre- 
sponding to de-coupled modes a^, has been given else- 
where [44j . and the generalization to include nonlineari- 
ties has been laid out in Ref. 0. Here we introduce cubic 
nonlinearities and make the rotating-wavc approximation 
(only terms with frequencies near wj, arc included in the 
equation of motion for afc) 0. This yields the following 
general coupled-mode equations: 
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As explained in Ref. 0, the nonlinear coefficients ay and 
/3fc depend on the specific geometry and materials, and 
express the strength of the nonlinear interactions. The 
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ajk terms describe self- and cross-phase modulation ef- 
fects which act to shift the cavity frequencies, while the 
f3k terms characterize the energy transfer (frequency con- 
version) between the modes. As noted in Ref. [U, these 
terms are constrained by energy conservation, which 
amounts to setting ^(|ao| 2 + \a m \ 2 + \a p \ 2 ) — (in the 
absence of external losses), yielding the following rela- 
tion: 

Wo/3q = LO m l3 m + UJ p f3 p . (5) 

In the following sections, for simplicity, we neglect losses 
such as linear absorption or radiation, i.e. we assume 
T s ,k = Tfe, and neglect nonlinear two-photon absorption, 



where is the electric field in the fcth mode and 
the denominators arise from the normalization of |<Zfc| 2 . 
As expected, Eqs. (JMZl) satisfy Eq. ((5), where Eq. © 
was obtained by imposing energy conservation on the 
TCMT equations without reference to the specific case 
of Maxwell's equations. 

There are six different ajk parameters [three SPM 
(etjj) and three XPM (ajk) coefficients], and in general 
they will all differ. However, from Eqs. (HH3) we see that 
they are all determined by similar modal integrals, lead 
to frequency shifting of the cavity frequencies, and all 
scale as 1/V, where V denotes a modal volume of the 
fields [4l|. Therefore, in the following sections, we begin 
by neglecting the frequency-shifting terms as in Ref. [l|, 
and then in Sec. IIV CI we study the essential effects of 
frequency shifting in the simplified case where all the co- 
efficients are equal (ajk = a). Of course, for a specific 
geometry one would calculate all coefficients Eqs. ([61 ITT)]) . 
but in this paper, we focus on the fundamental physics 
and phenomena rather than the precise behavior of a 



i.e. we assume aij are strictly real (two-photon absorp- 
tion effects can be minimized by selecting materials less 
susceptible to such processes). As noted below, these 
considerations do not qualitatively change our results, 
but merely act to slightly decrease the overall conversion 
efficiency once losses are included P, Q • 

The dependence of the coupling coefficients ay and 
fik on the geometry of the system can be obtained via a 
simple perturbative calculation involving the linear eigen- 
modes of the cavity, as described in Ref. 0- Carrying out 
this procedure to first order in \^ yields the following 
coupling coefficients: 
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specific geometry. 



III. QUANTUM-LIMITED VS. COMPLETE 
CONVERSION 

As described below, the DFWM process we consider 
here exhibits drastically different behavior depending on 
the ratio of Aw to ujq. In particular, there exist at 
least two distinct regimes of operation, corresponding to 
quantum-limited (|Aw| < wo) and complete (Ato > ljq) 
conversion. It turns out that, although our coupled-mode 
formalism is entirely classical, the same behaviors can 
be more easily understood by considering photon inter- 
actions in a quantum picture. Although this system is, of 
course, described by the general Manley-Rowe relations, 
which can be derived from both classical [3 UjJ [48[ and 
quantum [42|, [43[ arguments similar to those here, it is 
useful to review a basic picture of such limits and their 
physical consequences for the specific case of intra-cavity 
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Auj\ < Aui > u>o 




FIG. 2: (Color online) Diagram of nonlinear up-conversion 
process involving input light at u)o and uj m and output light at 
uip and uj m ■ The conversion efficiency of DFWM is determined 
by Auj, and photon energy conservation consideration (see 
text), leading to at least two different regimes of operation: 
(Left:) for |Aw| < ujo, two ujo pump photons and an signal 
uJm photon are converted into two uj m signal photons and 
an uj p photon. The input uj m photon is only necessary to 
initiate the conversion process and emerges unchanged after 
the interaction (indicted by red). (Right:) for Auj > ujo, two 
incoming ujo and a single uj m photon are combined to produce 
an cop photon. In contrast to the previous regime, the UJ m 
photon is energetically needed to produce the uj p photon. 

DFWM. 

Our focus in this manuscript is the up-conversion pro- 
cess (or interaction) corresponding to taking input light 
at frequencies ojq and oj m and generating output light at 
frequency uj p . Therefore, an appropriate figure of merit is 
the ratio of the output power in the oj p mode to the total 
input power, which we define as the absolute efficiency 

? ,= | SP ,_| 2 /(|S0, + | 2 + Kn, + | 2 ). 

As described in the previous section, the coupled- 
mode equations [Eqs. ([T]-|3])] follow from very general 
and purely classical considerations. The same consid- 
erations yield relationships between the frequencies and 
coupling coefficients of the problem, such as frequency 
conservation (uj m + uj p = 2ujq) and energy conservation 
(uJmPm + w p /3p = uiq/3q). Additional conservation rules 
that are perhaps best understood from quantum argu- 
ments, such as photon ener gy (?t|q;|) conservation and 
standard x^ selection rules [48|], also play a substantial 
role in the physics of nonlinear frequency conversion. In 
the case of the DFWM up-conversion process considered 
here, x' 3 ^ selection rules imply that nonlinear interac- 
tions can only be initiated if there exist at least three 
input photons: 2ujq photons and one uj m photon. 

In the | Aw | < ujq regime, there are at least two impor- 
tant features that can be understood from the above rela- 
tions: First, depletion of the signal input power (s„ lj+ ) is 
impossible, leading to a conversion efficiency rj < 1. Sec- 
ond, in order to maximize the total conversion efficiency, 
one desires s m> + to be as small as possible. These features 
can be understood by considering a simple picture of the 
nonlinear photon- photon interaction, as follows. From 
the DFWM x (3) selection rule [3, it follows that the 
creation of an oj p photon is accompanied by the destruc- 



tion of two ujq photons and one u> m photon. The latter, 
along with photon energy conservation, leads to the pro- 
cess considered in Fig. [5] (left), in which two ujq photons 
and an oj m photon interact to yield two uj m photons and 
an uj p photon. From the figure, and since 2ujq > uj p , one 
can see that the incident uj m photon (depicted in red) is 
merely required by the selection rule to initiate the 
interaction, and emerges unmodified, accompanied by an 
oj p photon and an additional uj m photon. Thus, it is 
clear that the input oj m photon does not actively partici- 
pate in the energy transfer and therefore merely reduces 
the maximum possible conversion efficiency. This implies 
that one desires a minimal input signal power to initiate 
the up-conversion. Effectively, the incident uj m photons 
are amplified by the conversion process (a similar ampli- 
fication effect is a crucial component in othe r nonline ar 
interactions, such as OPAs in x^ media [il . Il04 Il05l |). 
In addition, it is clear that complete depletion of the sig- 
nal photons, i.e. s mi _ = 0, is not possible for non-zero 
s m ,_|_, and therefore the conversion efficiency must be less 
than 100% (since the total input power is conserved). 
No such restriction is placed on s ,-, and therefore we 
expect that maximal efficiency will be obtained for ar- 
bitrarily low signal power and complete depletion of the 
pump power, i.e. So,- = 0. 

Based on these arguments, we can predict the max- 
imal efficiency of the conversion process by considering 
the ratio of the energy of the output oj p photon (hu) p ) 
to the energy of the three input photons [h(2ujQ + w m )]. 
Since the U) m photons can be provided with arbitrarily 
low amplitude, we therefore expect maximal efficiency to 
be achieved upon neglecting their contribution, i.e. we 
predict a maximal efficiency of: 

r hnax (\A U < Uo \) = ^ = ^-. (11) 

Note that this efficiency depends only on the ratio of Auj 
to ujo and h cancels, so it should appear in the classical 
limit as well. As we shall see in Sec. IIV Al this predic- 
tion is verified analytically by examining the steady-state 
solution of our coupled-mode equations. 

In the Auj > ujq regime, the conversion process is fun- 
damentally different and, in particular, complete deple- 
tion of the uj m and ujq photons is possible, leading to 
100% conversion efficiency. Basically, because uj p > 2ujq 
in this case, no additional photons are required to sat- 
isfy photon energy conservation, yielding the nonlinear 
interaction process depicted in Fig. [2] (right), where two 
input ujq photons and an uj m photon combine to pro- 
duce an uj p photon. Note that now the input uj m photon 
actively participates in the energy transfer, in contrast 
to the | Aw | < ujo regime, leading to a maximal conver- 
sion efficiency occurring when so,+ and s TOj + are both 
non-zero. Furthermore, since uj p is now the only prod- 
uct of the interaction, we expect that complete depletion 
of both the pump and signal powers, so._ = s mi _ = 0, 
should be possible, leading to 100% conversion efficiency. 
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As before, this can also be quantified by comparing the 
ratio of the output energy (hui p ) to the input energy 
[h(2uo + |w TO Q] (note that now the energy of w m pho- 
ton is h\oj m \), and the result follows from the fact that 
2w + \ui m \ = uj p . Again, we shall see in Sec. IIVB1 this 
prediction is validated analytically and directly from the 
couplcd-modc equations, yielding also the critical input 
powers at which 100% conversion is achieved. 

In this section, we made a number of predictions based 
on very general arguments relying on a quantum inter- 
pretation of the nonlinear interactions, allowing us to ob- 
tain predictions of maximal conversion efficiency. Our 
final results, of course, contained no factors of H and it 
is therefore not surprising that we recover the same re- 
sults (albeit with more detail, e.g. predictions of the 
values of critical powers) in the ensuing analysis of the 
purely-classical coupled-mode equations. Nevertheless, 
the heuristic quantum picture of Fig. [5] has the virtue of 
being simple and revealing, while the classical derivation 
is more complicated (although more quantitative). Sim- 
ilar quantum arguments have also proven useful in other 
context s, suc h as in many problems involving classical ra- 
diation |l06j |. or the recently-studied problem of op tical 
bonding/anti-bonding in waveguide structures [107j . 

IV. COUPLED-MODE ANALYSIS 

In order to gain a simple understanding of the system, 
we shall first consider frequency conversion in the ab- 
sence of self- and cross-phase modulation, i.e. atjk = 0. 
The nonzero-a case will be considered in Sec. IIV CI Sec- 
tion IIV Al focuses on the |Aw| < wo regime, whereas 
Sec. IIV Bl focuses on the Aw > wo regime. In both cases, 
we describe the solutions to the coupled-mode equations 
[Eqs. HH3]))] in the steady state, including the statbil- 
ity of these solutions and their dependence on the cavity 
parameters. 



A. |Aw| < Uo regime: Limited conversion 

Although the analysis in this section is general, for the 
purposes of plotting results we choose the specific pa- 
rameters: Ctjk = 0, To = T m = T p = 100/wo, P = 10~ 4 , 
and Aw = 0.05wo- The qualitative results remain un- 
changed as these parameters are varied, provided that 
the Q are large enough such that mode overlap is mini- 
mal as required by CMT. The influence of varying these 
parameters is discussed further at the end of the section. 

To understand the stability and dynamics of the non- 
linear coupled-mode equations in the quantum-limited 
regime, we apply the standard technique of identifying 
the fixed points of Eqs. ([T] [5]) and analyzing the sta bility 
of the linearized equations around each fixed point [108l | . 
A fixed point is given by a steady-state solution where 
the mode amplitudes vary as ak(t) = Ake ZUJkt , with the 
Ak being unknown constants. Plugging this steady-state 



ansatz into Eqs. (HH3|), we obtain three coupled polyno- 
mial equations in the parameters Aq, A m , A p , so,+j an d 
s m + . These polynomials were solved using Mathematica 
to obtain the mode energies |^4fc| 2 , which are then used 
to calculate the efficiency rj = |s p ,_| 2 /(|so.+ | 2 + |s m ,+ | 2 ). 
The phases of the A^ can be easily determined from the 
steady-state equations of motion; Aq and A m acquire the 
phases of so.+ and s m + respectively, while the phase of 
A p is that of fi v A\A* m rotated by 7r/2. Without loss of 
generality, so,+ an d s m> + can be chosen to be real. 

In general, this system has either one or three solu- 
tions, only one of which is ever stable. The stability 
and efficiency of this solution are shown in Fig. [3] for the 
specific parameters mentioned above. We observe that 
maximal conversion efficiency is obtained in the limit as 
input signal power s mi+ is reduced to zero, consistent 
with the discussion in the previous section. To obtain the 
maximum efficiency and the corresponding critical input 
powers, complete depletion of the pump (wo) photon is 
required, i.e. so.- =0 (note that one cannot require de- 
pletion of the signal photon, for the reasons discussed in 
the previous section). We find that the maximum effi- 
ciency 77, nax is obtained at Isf^+I 2 = -Po as |s m ,+ | 2 — > 0, 
where: 



To\/3o\\/T m T p \uJ m UJp\ 

top 1 / Aw\ 

W = ^ = 9 1 + > ( 13 ) 

2w 2 \ u>o J 

Note that Eq. (fT5|) is identical to the value predicted in 
the previous section. In the important case of narrow- 
band conversion, Aw| <C wo, the maximum efficiency 
is approximately 50%. (however, this is relative to the 
pump power — compared to the input signal alone, the 
output signal is amplified to an arbitrary degree). If Qo ~ 
Qm ~ Qp, then, as in THG [l|, the critical power scales 
as V/Q 2 , where V is the modal volume (recall that j3 ~ 
l/V). 

As Aw — > ujq, the maximum efficiency approaches 
unity, i.e. 100% conversion can be achieved in the limit. 
This limit is reminiscent of second-harmonic generation, 
since oj p = 2luq. However, the interaction process is fun- 
damentally different from the standard (x^) SHG in a 
number of ways. First, one is converting DC (w m w 0) 
light and wo pump light into 2wo- Second, the stabil- 
ity of this solution (described below) is quite different 
from that of SHG [iy, |47], [HJ . Finally, the critical power 
in this case, Po, diverges as l/yT — (Acj/cj ) 2 for Am 
near ujq. However, Aw close but not equal to ujq yields 
a reasonable P : for example, Alu = 0.95o;o yields effi- 
ciency rj = 0.975 with a critical power roughly three times 
the critical power for Aw near zero. Because this near- 
"SHG" situation involves coupling resonances at very dif- 
ferent frequency scales, it is reminiscent of using 
DFG to produce THz from infrared [60l ]. 

Eqs. (fT2HT3")) are only valid in the limit |s„ li+ | 2 — » 0, 
which is ideal from an efficiency perspective. How- 
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FIG. 3: (Color online) Color plot of the steady-state conver- 
se- | 2 /(|so,+ 1 2 + \s m ,+ \ 2 ) as a function of 
and |s mj + j 2 . for a system consisting of 



sion efficiency r\ 
input power \so,- 



Acj = 0.05ojo, P = 10 , and to = r p 



100. 



Both pow- 

|2 



ers are normalized by the critical power P c (|s mj +| — > 0) = 
Po = 2/to\P\\/ T m T p \oj m u)p | (black dot). The shaded region in- 
dicates that the solution is unstable. The curves P± indicate 
the powers at which depletion of the ujo input light is achieved, 
i.e. so,- = 0; the critical power P c (|s m ,+ | 2 ) is defined as the 
total input power that yields the highest stable efficiency for 
any given |s m> +| 2 . The dash line is the cross-section shown in 
Fig- HI 
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FIG. 4: (Color online) Bifurcation diagram of the steady-state 
efficiency r], normalized by the quantum-limited maximum ef- 



ficiency r/ n 



o p /loo, 



as a function of 

2 



|so,+ | , normalized 

by Pq, for signal power |s m ,+ | 2 = O.lPo (indicated by the 
black dashed line of Fig. [3]). Red/blue correspond to a sta- 
ble/unstable solution (note that the two bifurcating solutions 
are always unstable). The green dashed line illustrates the 
bounds of the limit cycles obtained from time domain simu- 
lations, where the solid green line yields the average over the 
cycle. (Inset:) Efficiency as a function of time in units of 
the period T p = 2-k/uj p in a regime where there exists a limit 
cycle. 



ever, it is interesting to consider the system for non- 
infinitesimal s mi _|_, in which case we solve for the input 
power that yields a stable solution with maximal effi- 
ciency for a given s m + . We denote this input power by 
Pc{\s m ,+ \ 2 ) = K+\ d + Iwl 2 , where |s^+| 2 (a func- 
tion of |s mj _)-| 2 ) is defined to be the pump power required 
to achieve maximum, stable conversion efficiency for a 
given signal power |s m , + | 2 . As seen in Fig. [3l this effi- 
ciency is always < J7 m ax, and P c Pq as s m ,+ — > 0. In 
the non-zero |s m + | 2 regime, P c does not correspond to 
complete depletion of the pump. Requiring pump deple- 
tion (so,- = 0) for a given signal power |s mi +| 2 yields two 
pump powers, which we label P±{\s m .+ \ 2 )- P+ (|s m .+ | 2 ) 
does indeed provide a solution with maximal efficiency, 
however this solution is always unstable. As seen from 
Fig. |31 only for small signal power s m + does depletion of 
the pump lead to maximal efficiency. 

In general, to obtain the largest efficiency while retain- 
ing stability, one would aim to operate with low signal 
power |s m + | 2 and use a pump power near the critical 
power Po given in Eq. (|12[) . However, it is interesting to 
consider the unstable solutions, because they turn out to 
be related to limit cycles. As mentioned above, the sys- 
tem contains either one or three steady-state solutions 
for given input powers. Fig. [4] plots these stable and un- 
stable solutions as a function of pump power |s ,+ | 2 at 
fixed signal power |s m + | 2 = O.lPo, corresponding to the 
horizontal dashed line in Fig. [31 For low input pump 



power |so.+ | 2 , the system has a single steady-state solu- 
tion; as the pump power is increased, the system expe- 
riences a bifurcation yielding two unstable solutions. As 
mentioned above, the higher efficiency solution emerging 
from the bifurcation achieves a maximum corresponding 
at |so,+ | 2 = P+, coinciding with complete depletion of 
the pump (so,- = 0), but this maximal efficiency solu- 
tion is always unstable; note that there may be a stable 
solution at |so.+ | 2 = P+, but the stable solution will have 
a lower efficiency than the maximal, unstable solution, as 
shown in Fig. 01 Furthermore, the original stable solu- 
tion eventually becomes unstable as the pump power is 
increased (this can occur before or after the bifurcation, 
depending on the system parameters); this onset of in- 
stability coincides with the onset of limit cycles, stable 
oscillating-cfficicncy solutions. An example of these limit 
cycles are shown in Fig. 21 where the green dashed lines 
indicate the bounds of the oscillations and the solid green 
line gives the average. The limit cycles are plotted as a 
function of time in the inset of Fig. 2) The limit cycles 
shown here were obtained by numerically time-evolving 
the coupled-mode equations. In general, we find that 
these limit cycles oscillate with a period proportional to 

Tp. 

Figures [3HH describe a system corresponding to a par- 
ticular set of values for the parameters Aw and . Qual- 
itatively, the most important features of the figures re- 
main largely unchanged as these parameters are var- 



8 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.! 

\s m ,+ \ 2 /Po 

FIG. 5: (Color online) Plot of the steady-state efficiency 
77 (solid lines) along the critical solution [total input power 
Pc{\s m ,+ \ 2 ) = s m ,+ | 2 + I s§, -1- 1 2 that yields the maximum ef- 
ficiency for a given s m ,+ | 2 (solid white curve of Fig. E}] and 
the value of P c (dashed lines) as a function of s mj +| 2 , nor- 
malized by Po, for three different values of Aoj: O.lwo (red), 
0.5wo (blue), and 0.9u;o (green). The kinks in the P c curves 
correspond to the point U where P c reaches the region of in- 
stability (see Fig. |3]) . The coupling lifetimes r and coefficient 
/3 of the system are equivalent to those of Fig. [3] 



ied. Basically, there exist at most three solutions to the 
coupled-mode equations, one of which has a finite region 
of stability as a function of ,so,+ and s m!+ , with the gen- 
eral shape that is shown in Fig. |3l and two others that 
are always unstable and bifurcate at a finite so,+- There 
are however, some differences to note: First, as Aw in- 
creases from 0, the maximum steady-state efficiency also 
increases, asymptoting to rj = 1 as Aw — > loq. This was 
obtained analytically and is quantified in Eq. (|13p . Un- 
fortunately, we find that as Aw increases, the region of 
instability in Fig. [3] also increases, and furthermore, the 
conversion efficiency at finite s lri!+ also drops off more 
rapidly. (In particular, we observe in the "SHG" limit of 
Aw — > wo, the system becomes largely unstable except for 
very low signal powers.) These tendencies arc depicted 
in Fig. [SJ which plots P c (|s mj +| 2 ) and the corresponding 
conversion efficiency for different values of Aw. The kinks 
observed in the plots of P c are due to the discontinuity 
in the slope of the P c curve as it reaches the region of 
instability, corresponding to the point U in Fig. [3] 

Varying Tk does not affect the maximum possible ef- 
ficiency and also leaves Fig. [3] qualitatively unchanged, 
changing only the scale of the critical input power Pq. 
The stability of the system however, does depend on the 
relative lifetimes of the cavity modes. In particular, the 
stability depends largely on the ratio to/t p , and decreases 
weakly as r rn increases with respect to either To or r p . 
This makes sense since, as argued in Sec. IIII1 the w m 
photons do not actively participate in the energy trans- 
fer. (A similar dependence on the ratio of the lifetimes 



was also observed in the case of THG [l|.) More quanti- 
tatively, we follow the position of the point U (the point 
where P c reaches the region of instability) as the are 
varied. Assuming equal modal lifetimes (tq = r m = T p 
as in Fig. [3]), we find that U lies at critical input powers 
,s ,+ | 2 ~ 1.28Po an d |s m .+ | 2 « 0.35p). Increasing t /t p , 
from 1 to 10, we find that U moves to so,+ | 2 ~ 10Po and 
w 4.75P . However, if we instead keep r = r p 
and increase r m such that T m /ro = T m /r p = 10, U moves 
only to |s ,+ | 2 « 1.05P and s mj+ | 2 w 0.27P . Note 
that, as mentioned previously, maximal stable conver- 
sion efficiency is obtained for low signal power |s TOj +| 2 
and input power |so.+ | 2 near the critical power Po, re- 
gardless of Tfc. We note that rescaling /3 simply scales 
the input power and therefore changing f3 does not affect 
the dynamics. 

Thus far, we have focused on the up-convcrsion pro- 
cess: taking input light at frequencies wo and w m and gen- 
erating output light at frequency w p > wo. However, it 
suffices to consider the above system when Aw < to un- 
derstand the physics of the alternative, down-conversion 
process: taking input light at frequencies wo and w p and 
generating output light at frequency w m . For Aw < 0, 
we effectively have w m -H- w p . In this regime, all of the 
above analysis holds, and in particular, the maximal effi- 
ciency, given by Eq. (fT3"]) . is obtained as |s m ,+ | 2 — > with 
I so. + 1 2 = Po- Similarly, the stability of the solutions fol- 
low similar trends to those outlined above. 



B. Auj > ixjq regime: Complete conversion 



When Aw is larger than wo, we argued in Sec. Illll that 
the system is capable of complete conversion, i.e. rj = 1. 
In this section, we demonstrate the existence of a critical 
steady-state solution to the classical coupled mode equa- 
tions with complete conversion and analyze the stability 
of this critical solution, as well as relate DFWM to our 
previous work on THG [H, 0] • 

As in the previous section, we consider the equations 
of motion Eqs. ([TH3]) in the steady state. To obtain the 
critical solution, we again require depletion of the pump 
power, i.e. so.- = 0. However, as argued in Sec. lIIIi com- 
plete depletion of the signal, s mi _ = must also occur. 
Recall from Sec. Illll that complete w m depletion is possi- 
ble in the Aw > wo regime since the up-conversion pro- 
cess does not produce w m photons (see Fig. [5]). Imposing 
the depletion constraints on the steady-state equations of 
motion yields the following critical cavity energies |a"' lt: | 2 : 



(14) 

(15) 
(16) 



crit 


2 = 


1 




\Pm \y/T m T p W m W p | 


a crit 
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T ™l W m l„crit|2 






2r wo |a ° 1 ' 
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:■ _ 
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FIG. 6: (Color online) Plot of the critical powers |so,+ | 2 
(blue), |s ?7li +| 2 (red), and maximum steady-state efficiency 
r\ (green) as a function of Aw/wo (the tilde over the criti- 
cal powers indicates that the values have been rescaled by 
the factor 4/ro|/3o| v^W^Wo). The vertical dashed lines at 
Aw = wo and Aw = 2wo indicate special degenerate regimes, 
corresponding to "second harmonic generation" (SHG) and 
third harmonic generation (THG). (Note the discontinuity in 
|s m ,+ | 2 located at Aw — 2wq, explained in the text). 

which lead to the following critical powers: 

K+l 2 = (17) 
K+| 2 = ^o, (18) 

where Pq is given by Eq. (fT^j) . Solving for the correspond- 
ing output signal s p ._| 2 , the output power is indeed 
100% of the input power, as required by energy conserva- 
tion. (In contrast, the assumption that so.- = s m .- = 
in the Aw| < ujq case yields no solution). Note that 
the critical signal power |s™+| 2 is now non-zero, due to 
the fact that the energy from the signal uj m photons is 
necessary to produce the output w p photons. This is in 
contrast with the |Aw| < ujq regime where maximal con- 
version efficiency was only achieved in the limit as input 
signal power |s TOj+ | 2 decreased to zero. The critical pump 
and signal powers, with the corresponding maximum ef- 
ficiency ?7, are plotted versus Aw in Fig. [5] for both Aw 
regimes. 

As may be noted from Fig. |5J there are two partic- 
ular values of Aw that warrant special attention when 
Alu > ujq. The first case, when Auj = ljq, the "SHG" 
case, was discussed in the previous section. The sec- 
ond case is when Alo = 2loq. In this case, u> m = —luq 
and LOp = 3wo, reminiscent of third- harmonic generation 
(THG). In fact, this case of DFWM corresponds exactly 
to x THG, and thus Aw > wo strictly generalizes our 
previous THG analysis [l|. To see this, some care must 
be taken to adjust the coupling coefficients 8k given in 
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FIG. 7: (Color online) Stability contours (number of stable 
solutions) as a function of modal lifetimes r m and r p , normal- 
ized by to, pumping at the critical input powers |so'+| 2 and 
|Sm'+| 2 . The stability in the Aw > wo regime is independent 
of the value of Aw. 



Eqs. ([6]-[7]) to properly implement the rotating wave ap- 
proximation; since u> m = — w , we have that a m = clq, 
and thus B -> B + B* m and B m B m + Bq . This results 
in fio = B* n = 3/3*, exactly as shown in [l(. Furthermore, 
we have |sg^| 2 = \s c ^ + \ 2 = Po (note that this differs 
by a factor of two from Eq. (fT8|) . due to the adjusted /3k 
values); upon requiring that to = r m , this recovers the 
critical power previously obtained for THG [l[ . Note that 
the correspondence between Aw = 2wo and THG is 
exact, whereas the Aw = wo limit has little in common 
with x^ 2 ) SHG as discussed above. 

The existence of an so,- = s mj _ = solution having 
demonstrated the existence of critical powers where 100% 
conversion can be achieved, we are now interested in char- 
acterizing the system at this critical power by studying 
all of the fixed points. These fixed points were obtained 
using Mathematica as in the previous section, and their 
stability was determined via linear stability analysis as 
before. For the critical input power, the steady-state 
equations of motion yield three solutions; however, in 
contrast to the |Aw| < ujq regime, there exists multi- 
stability when Aw > wo. Similar to the case of THG 
(Aw = wo), the system is either singly stable, doubly 
stable, or unstable, depending on the values of the mode 
lifetimes (see Fig. [7]). In this Aw > wo regime, the 
stability of the solutions does not depend on Aw, again 
in contrast with the quantum-limited regime. Unlike the 
| Aw | < wo regime, the value of r m now plays a significant 
role in the stability of the solutions. 
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C. Self- and Cross-Phase Modulation (a ^ 0) 



Finally, we briefly consider the effects of SPM and 
XPM. This corresponds to taking the coefficients ajk to 
be non-zero; as mentioned above, for simplicity, we take 
all the coefficients to be equal, i.e. ajk = a for all j, k. 
The main effect of SPM and XPM in Eqs. (TH^ is to 
shift the resonant frequencies of the cavity in proportion 
to the energy of the modes in the cavity. Generally, this 
drives the frequency input light off resonance and there- 
fore degrades the overall conversion efficiency obtained 
in Sees. JTVAl and IIVB1 as shown in Ref. [l|. However, 
in Ref. [l|, we showed that one simple way to overcome 
this difficulty is to pre-shift the cavity resonant frequen- 
cies so as to compensate for the SPM/XPM effects when 
operating near the critical input power. Unfortunately, 
this will inevitably affect the stability analysis obtained 
in the a = case, and therefore a new analysis that 
includes SPM/XPM effects must be performed. In this 
remainder of this section, we only analyze the stability 
of the maximal- efficiency solutions obtained in Secs. llV Al 
and IIV B\ and in particular, we find that 100% photon- 
conversion efficiency can be obtained in this case as well. 



The change in cavity frequency due to SPM/XPM can 
be accounted for by a pre-shifting technique described in 
Ref.[l|. In particular, the a terms in Eqs. (PES) act to shift 
the cavity resonant frequencies from w£ av — > L , spoil- 
ing the frequency-conservation relations necessary for ef- 
ficient nonlinear frequency conversion as well as detuning 
the resonances from the input light. However, one can 
simply design the cavity frequencies to be resonant at the 



shifted frequencies, i.e. w£ av = w£ , for a given steady- 
state solution. For the critical solutions corresponding to 
100% photon-conversion efficiency, this implies that the 
new cavity frequencies will be given by [l[ : 
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where la™ 1 '! 2 are the energies of the modes at critical 
power. For cavities resonances w£ av , the new equations 



of motion are given by: 

iw§ av (l - a o|a | 2 - a 0m \a m \ 2 - ao P \a p \ 2 ) 
a - iu} /3 a aQa m a p 



~~dt 
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iw!; av (l - a p0 \a \ - a pm \a m \ ~ a pp \a p \ ) 



a p - iui p j3 p ala* m , (24) 



Note that the frequencies u>k multiplying the terms 
do not need to be shifted, since the terms introduced by 
such a shifting will be higher order in x*- 3 -*. By inspec- 
tion, we observe that the solutions obtained in Sees. IIV Al 
and IIV Bl at critical input power a£ rlt are also solutions 
of Eqs. (|2"2Tl2"5)) . but as explained above, their stability 
may change. Using the results from Sees. IIV XI and HVB| 
we now study the stability properties of these solutions 
in the two Aw regimes. 

We first consider the Alu < ljq regime. As in Sec. IIV Al 
we restrict our analysis to a specific parameter regime 
( To = r m = t p = 100/wo, /3 = 10" 4 , and Auj = 0.05w ) 
for simplicity, although our qualitative conclusions ap- 
ply to other parameter ranges. As discussed above in 
Sec. IIV A| the maximal efficiency is obtained for input 
light with |so,+ | 2 = Po as |s mi+ | 2 — > 0. Since one must 
always pump with finite |s mj _|_| 2 , and there are no ana- 
lytic solutions in this case, we solve for the field energies 
l a fc rlt | 2 numerically at a small |s m .+ | 2 and for |so,+ | 2 = Po 
in the case of a = in order to compute the shifted 
frequencies Eqs. P^H2T|) . This allows us to solve the 
coupled-mode equations Eqs. (|22ll23)) and therefore ob- 
tain the steady-state field amplitudes and phases. As in 
Ref. Q], the inclusion of self- and cross-phase modulation 
introduces new steady-state solutions absent in the a = 
case, and the stability of the old and new solutions are 
then examined again via a linear stability analysis, as 
in Sec. IIV Al In particular, we find that the inclusion of 
SPM/XPM does not destroy the stability of the maximal 
efficiency solution in the a = case studied in Sec. IIV Al 
and in fact creates additional stable solutions, as shown 
in Fig. E 

A similar analysis can be performed in the Aw > 
wo regime, where it is possible to obtain the analytic 
form of the maximal efficiency solutions [Eqs. (fT^l fTB|) in 
Sec. IIVB) . Wc find that, as in the previous regime, the 
presence of a introduces additional stable solutions, while 
retaining the original 100% efficiency a = solution, over 
finite regions of the parameter space. 

The presence of SPM/XPM in our system provides an 
opportunity to observe rich and interesting dynamical 
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FIG. 8: (Color online) Contour plot of number of stable so- 
lutions (n a ) as a function of a//3 and jso,+ | 2 /Pu, for input 
pump power |s mj+ | 2 = O.lPo, and for the system described in 
the text. 

behaviors, including limit cycles and hysteresis effects, 
that we do not explore in this paper. As noted in this 
section, the inclusion of these effects is not prohibitive for 
100% nonlinear frequency conversion although predicting 
which parameter regimes allow for such conversion will 
depend on the system under question. In the future, we 
plan to examine SPM/XPM effects in more detail for 
realistic geometries with realistic values of ay and /3j. 
As in Ref. fl the presence of multiple stable solutions 
means that the manner in which the source is initiated 
will determine which solution is excited, but a simple 
initialization procedure similar to that in Ref. [l] should 
be possible to excite the maximal-efficiency solution. 

V. CONCLUSION 

By exploiting a simple but rigorous couplcd-mode the- 
ory framework, we have demonstrated the possibility of 



achieving highly-efficient (low-power) DWFM in triply- 
resonant cavities, similar to our previous work in SHG 
and THG 0, Q . We conclude that there are two main 
regimes of operation, determined by the ratio of Aw to 
ujq. In particular, whereas the maximal efficiency obtain- 
able in the Alu < ujq regime, corresponding to conversion 
between closely-spaced resonances, is bounded above by 
a quantum-limited process, there is no such bound when 
Auj > loq. In both regimes, a suitable choice of system 
parameters leads to stable, maximal-efficiency nonlinear 
frequency conversion, even in the presence of SPM and 
XPM effects. We remark that all of the results obtained 
in this paper correspond to the idealized case of loss- 
less interactions, since the main focus of the paper is in 
examining the basic considerations involved in operat- 
ing with these systems rather than predicting results for 
specific experimentally-relevant systems. Nevertheless, 
based on our previous experience with SHG and THG 
P, 0, we expect that linear and nonlinear losses, e.g. 
coming from radiation or material absorption, will only 
act to slightly decrease the overall conversion efficiency 
and will not affect the qualitative predictions here. In a 
future manuscript, we plan to explore DFWM in a re- 
alistic geometry such as a ring resonator coupled to an 
index-guided waveguide and study some of the dynamical 
effects arising from SPM/XPM. 



Acknowledgments 



This work was supported by the MRSEC Program 
of the NSF under Award No. DMR-0819762 and by 
the U.S.A. R.O. through the ISN under Contract No. 
W911NF-07-D-0004. We are also grateful to Jorge 
Bravo- Abad at U. Autonoma de Madrid for helpful dis- 
cussions. 



[1] H. Hashemi, A. W. Rodriguez, J. D. Joannopoulos, 
M. Soljacic, and S. G. Johnson, Phys. Rev. A 79, 013812 
(2009). 

[2] A. Rodriguez, M. Soljacic, J. D. Joannopulos, and S. G. 
Johnson, Opt. Express 15, 7303 (2007). 

[3] B. E. A. Saleh and M. C. Teich, Fundamentals of Pho- 
tonics (Wiley, 1991). 

[4] J. U. Fiirst, D. V. Strekalov, D. Elser, M. Lassen, U. L. 
Andersen, C. Marquardt, and G. Leuchs, Phys. Rev. 



Lett. 104, 153901 (2010). 
[5] M. Bieler, IEEE J. Sel. Top. Quant. Elec. 14, 458 
(2008). 

[6] R. E. Hamam, M. Ibanescu, E. J. Reed, P. Bermel, S. G. 
Johnson, E. Ippen, J. D. Joannopoulos, and M. Soljacic, 
Opt. Express 12, 2102 (2008). 

[7] I. B. Burgess, Y. Zhang, M. W. McCutcheon, A. W. Ro- 
driguez, J. Bravo-Abad, S. G. Johnson, and M. Loncar, 
Opt. Express 17, 20099 (2009). 



12 



[8] J. Bravo-Abad, A. W. Rodriguez, J. D. Joannopoulos, 
P. T. Rakich, S. G. Johnson, and M. Soljacic, Appl. 
Phys. Lett. 96, 101110 (2010). 
[9] J. Hald, Optics Communications 197, 169 (2001). 
[10] R. Lifshitz, A. Arie, and A. Bahabad, Phys. Rev. Lett. 

95, 133901 (2005). 
[11] B. Freedman, G. Bartal, M. Segev, R. Lifshitz, D. N. 
Christodoulides, and J. W. Fleischer, Nature 440, 1166 

(2006) . 

[12] Y. A. Morozov, I. S. Nefedov, V. Y. Aleshkin, and I. V. 

Krasnikova, Semiconductors 39, 113 (2005). 
[13] T. J. Kippenberg, S. M. Spillane, and K. J. Vahala, 

Phys. Rev. Lett. 93, 083904 (2004). 
[14] P. Bermel, A. Rodriguez, J. D. Joannopoulos, and 

M. Soljacic, Phys. Rev. Lett. 99, 053601 (2007). 
[15] J. Bravo-Abad, A. W. Rodriguez, P. Bermel, S. G. John- 
son, J. D. Joannopoulos, and M. Soljacic, Opt. Express 

15, 16161 (2007). 
[16] M. M. Fejer, Phys. Today 47, 25 (1994). 
[17] V. S. Ilchenko, A. A. Savchenkov, A. B. Matsko, and 

L. Maleki, Phys. Rev. Lett. 92, 043903 (2004). 
[18] F. S. Felber and J. H. Marburger, Appl. Phys. Lett. 28, 

731 (1976). 

[19] A. Parini, G. Bellanca, S. Trillo, M. Conforti, A. Lo- 
catelli, and C. De Angelis, J. Opt. Soc. Am. B 24, 2229 

(2007) . 

[20] H. M. Gibbs, Optical Bistability: Controlling Light with 
Light (Academic Press, Orlando, FL, 1985). 

[21] A. R. Cowan and J. F. Young, Phys. Rev. E 68, 046606 
(2003). 

[22] M. Soljacioc, M. Imabenscu, S. G. Johnson, J. D. 

Joannopoulos, and Y. Fink, Opt. Lett. 28, 516 (2003). 
[23] M. F. Yanik, S. Fan, and M. Soljacic, Appl. Phys. Lett. 

83, 2739 (2003). 
[24] G. S. Dutta and J. Jolly, Pranama-journal of Physics 

50, 239 (1988). 
[25] Q. F. Xu and M. Lipson, Opt. Lett. 31, 341 (2006). 
[26] M. Soljacic, M. Ibanescu, S. G. Johnson, Y. Fink, and 

J. D. Joannopoulos, Phys. Rev. E Rapid Commun. 66, 

055601(R) (2002). 
[27] E. Centeno and D. Felbacq, Phys. Rev. B 62, R7683 

(2000). 

[28] T. Tanabe, M. Notomi, S. Mitsugi, A. Shinya, and 
E. Kuramochi, Opt. Lett. 30, 2575 (2005). 

[29] M. Notomi, A. Shinya, S. Mitsugi, G. Kira, E. Ku- 
ramochi, and T. Tanabe, Opt. Express 13, 2678 (2005). 

[30] A. Dorsel, J. McCullen, P. Meystre, E. Vignes, and 
H. Walther, Phys. Rev. Lett. 51, 1550 (1983). 

[31] K. Y. Billah and M. Shinozuka, Phys. Rev. A 42, 7492 
(1990). 

[32] K. J. Vahala, Nature 424, 839 (2003). 

[33] M. Ferrera, D. Duchesne, L. Razzari, M. Peccianti, 

R. Morandotti, P. Cheben, S. Janz, D. X. Xu, B. E. 

Little, S. Chu, et al., OE 17, 14098 (2009). 
[34] M. Ferrera, L. . Razzari, D. Duchesne, R. Morandotti, 

Z. Yang, M. Liscidini, J. E. Sipe, S. Chu, B. E. Little, 

and D. J. Moss, Nature Photonics 2, 737 (2008). 
[35] A. C. Turner, M. A. Foster, A. L. Gaeta, and M. Lipson, 

OE 16, 4881 (2008). 
[36] J. E. Heebner, N. N. Lepeshkin, A. Schweinsberg, G. W. 

Wicks, R. W. Boyd, R. Grover, and P. T. Ho, OL 29, 

769 (2004). 

[37] P. P. Absil, J. V. Hryniewicz, B. E. Little, P. S. Cho, 
R. A. Wilson, L. G. Joneckis, and P. T. Ho, OL 25, 554 



(2000). 

[38] I. H. Agha, Y. Okawachi, M. A. Foster, J. E. Sharping, 
and A. L. Gaeta, PRA 76, 043837 (2007). 

[39] D. H. Broaddus, M. A. Foster, I. H. Agha, J. T. Robin- 
son, M. Lipson, and A. L. Gaeta, Opt. Express 17, 5998 

(2009) . 

[40] D. Duchesne, L. Razzari, M. Ferrera, R. Morandotti, 
S. Chu, B. Little, and D. Moss, in Optical Fiber Commu- 
nication (OFC), collocated National Fiber Optic Engi- 
neers Conference, 2010 Conference on (OFC/NFOEC) 

(2010) , pp. 1 -3. 

[41] J. D. Joannopoulos, S. G. Johnson, J. N. Winn, and 
R. D. Meade, Photonic Crystals: Molding the Flow of 
Light (Princeton University Press, 2008), 2nd ed., URL 
http : //ab-initio .mit . edu/book 

[42] M. T. Weiss, Proc. IRE 45, 1012 (1957). 

[43] J. Brown, EL 1, 23 (1965). 

[44] H. A. Haus, Waves and Fields in Optoelectronics 

(Prentice-Hall, Englewood Cliffs, NJ, 1984). 
[45] H. A. Haus and W. Huang, Proc. of the IEEE 79, 1505 

(1991) . 

[46] P. D. Drummond, K. J. McNeil, and D. F. Walls, Optica 

Acta. 27, 321 (1980). 
[47] K. Grygiel and P. Szlatchetka, Opt. Comm. 91, 241 

(1992) . 

[48] R. W. Boyd, Nonlinear Optics (Academic Press, Cali- 
fornia, 1992). 

[49] H. M. Gibbs, G. Khitrova, and N. Peyghambarian, Non- 
linear Photonics (Spring- Verlag, 1990). 

[50] R. G. Smith, IEEE J. Quantum Electron. 6, 215 (1970). 

[51] L.-A. Wu, M. Xiao, and H. J. Kimble, JOSA-B 4, 1465 
(1987). 

[52] Z. Y. Ou and H. J. Kimble, Opt. Lett. 18, 1053 (1993). 

[53] R. Paschotta, K. Fiedler, P. Kurz, and J. Mlynek, Appl. 
Phys. Lett. 58, 117 (1994). 

[54] V. Berger, J. Opt. Soc. Am. B 14, 1351 (1997). 

[55] I. I. Zootoverkh, K. N. V., and E. G. Lariontsev, Quan- 
tum Electronics 30, 565 (2000). 

[56] B. Maes, P. Bienstman, and R. Baets, J. Opt. Soc. 
Am. B 22, 1378 (2005). 

[57] M. Liscidini and L. A. Andreani, Phys. Rev. E 73, 
016613 (2006). 

[58] S. Schiller, Ph.D. thesis, Stanford University, Stanford, 
CA (1993). 

[59] A. Ashkin, G. D. Boyd, and J. M. Dziedzic, IEEE J. 

Quantum Electron. 2, 109 (1966). 
[60] I. B. Burgess, A. W. Rodriguez, M. W. McCutecheon, 

J. Bravo-Abad, Y. Zhang, S. G. Johnson, and M. Lon- 

car, OE 17, 9241 (2009). 
[61] M. W. McCutcheon, D. E. Chang, Y. Zhang, M. D. 

Lukin, and M. Loncar, arXiv:0903.4706 (2010). 
[62] Y. H. Avetisyan, Proc. SPIE 3795, 501 (1999). 
[63] R. A. Baumgartner and R. L. Byer, IEEE J. Quantum 

Electron. 15, 432 (1979). 
[64] A. Yariv, Quantum Electronics, 3rd ed. (Wiley, New 

York, 1998). 

[65] G. T. Moore, K. Koch, and E. C. Cheung, Optics Com- 
munications 113, 463 (1995). 

[66] M. Liscidini and L. A. Andreani, Appl. Phys. Lett. 85, 
1883 (2004). 

[67] Y. Dumeige and P. Feron, Phys. Rev. A 74, 063804 
(2006). 

[68] S. Stivala, A. C. Busacca, A. Pasquazi, R. L. Oliveri, 
R. Morandotti, and G. Assanto, OL 35, 363 (2010). 



13 



A. I. Ferguson and M. H. Dunn, IEEE J. Quantum Elec- 
tron. 13, 751 (1977). 

M. Brieger, H. Busener, A. Hese, F. V. Moers, and 
A. Renn, Opt. Commun. 38, 423 (1981). 
J. C. Bergquist, H. Hemmati, and W. M. Itano, Opt. 
Commun. 43, 437 (1982). 

W. J. Kozlovsky, W. P. Risk, W. Lenth, B. G. Kim, 
G. L. Bona, H. Jaeckel, and D. J. Webb, Appl. Phys. 
Lett. 65, 525 (1994). 

G. J. Dixon, C. E. Tanner, and C. E. Wieman, Opt. 
Lett. 14, 731 (1989). 

M. J. Collet and R. B. Levien, Phys. Rev. A 43, 5068 
(1991). 

M. A. Persaud, J. M. Tolchard, and A. I. Ferguson, 
IEEE J. Quantum Electron. 26, 1253 (1990). 
K. Schneider, S. Schiller, and J. Mlynek, Opt. Lett. 21, 
1999 (1996). 

X. Mu, Y. J. Ding, H. Yang, and G. J. Salamo, Appl. 
Phys. Lett. 79, 569 (2001). 

G. McConnell, A. I. Ferguson, and N. Langford, J. Phys. 
D: Appl.Phys 34, 2408 (2001). 

T. V. Dolgova, A. I. Maidykovski, M. G. Martemyanov, 
A. A. Fedyanin, O. A. Aktsipetrov, G. Marowsky, V. A. 
Yakovlev, G. Mattei, N. Ohta, and S. Nakabayashi, 
J. Opt. Soc. Am. B 19, 2129 (2002). 
T.-M. Liu, C.-T. Yu, and C.-K. Sun, Appl. Phys. Lett. 
86, 061112 (2005). 

L. Scaccabarozzi, M. M. Fejer, Y. Huo, S. Fan, X. Yu, 
and J. S. Harris, OL 31, 3626 (2006). 

C. M. Savage and D. F. Walls, Optica Acta 30, 557 
(1983). 

E. Abraham, W. J. Firth, and J. Carr, Physics Lett. A 
pp. 47-51 (1982). 

A. Siegman, Lasers (University Science Books, Mill Val- 
ley, CA, 1986). 

M. D. Levenson, R. M. Shelby, A. Aspect, M. Reid, and 

D. F. Walls, Phys. Rev. A 32, 1550 (1985). 

N. Shibata, R. P. Braun, and R. G. Waarts, IEEE J. 

Quantum Electron. QE-23, 1205 (1987). 

K. Inque, J. Lightwave Tech. 10, 1553 (1992). 

M. Karisson, JOSA-B 15, 2269 (1998). 

J. Hansryd, P. A. Andrekson, M. Westlund, J. Li, and 

P. O. Hedekvist, IEEE J. Quantum Electron. 8, 506 



(2002). 

[90] L. Deng, E. W. Hagley, M. Trippenbach, Y. Band, P. S. 
Julienne, J. E. Simsarian, K. Helmerson, S. L. Rolston, 
and W. D. Phillips, Nature 398, 218 (1999). 

[91] J. Feinberg, Opt. Lett. 7, 486 (1982). 

[92] S. Watanabe, T. Naito, and T. Chikama, PTL 5, 92 
(1993). 

[93] G. P. Agrawal, Fiber-Optic Communication Systems 

(Wiley, Hoboken, NJ, 2002). 
[94] H. P. Yuen and J. H. Shapiro, Opt. Lett. 4, 334 (1979). 
[95] R. E. Slusher, L. W. Hollberg, B. Yurke, J. C. Mertz, 

and J. F. Valley, Phys. Rev. Lett. 55, 2409 (1985). 
[96] A. Pasquazi, Y. Park, J. A. na, F. Legare, R. Moran- 

dotti, B. E. Little, S. T. Chu, and D. J. Moss, Opt. 

Express 18, 7634 (2010). 
[97] A. Pasquazi, R. Ahmad, M. Rochette, M. Lamont, B. E. 

Little, S. T. Chu, R. Morandotti, and D. J. Moss, Opt. 

Express 18, 3858 (2010). 
[98] G. Bartal, O. Manela, and M. Segev, Phys. Rev. Lett. 

97, 073906 (2006). 
[99] Y. C. Huang and Y. Y. Lin, JOSA-B 21, 777 (2003). 
[100] M. Troccoli, A. Belyanin, F. Capasso, E. Cubukcu, D. L. 

Sivco, and A. Y. Cho, Nature 433, 845 (2005). 
[101] S. Maas, Nonlinear microwave and RF circuits (Artech 

House, 2003). 

[102] A. Taflove and S. C. Hagness, Computational Electro- 
dynamics: The Finite- Difference Time-Domain Method 
(Artech, Norwood, MA, 2000). 

[103] W. Suh, Z. Wang, and S. Fan, IEEE J. Quantum Elec- 
tron. 40, 1511 (2004). 

[104] J. A. Armstrong, N. Bloembergen, J. Ducuing, and P. S. 
Pershan, Phys. Rev. 127, 1918 (1962). 

[105] S. K. Choi, R. D. Li, C. Kim, and P. Kumar, JOSA-B 
14, 1564 (1997). 

[106] J. D. Jackson, Classical Electrodynamics (Wiley, New 
York, 1998), 3rd ed. 

[107] M. L. Povinelli, M. Loncar, M. Ibanescu, E. J. Smythe, 
S. G. Johnson, F. Capasso, and J. D. Joannopoulos, 
Opt. Lett. 30, 3042 (2005). 

[108] M. Tabor, Chaos and Integrability in Nonlinear Dynam- 
ics: An Introduction (Wiley, New York, 1989). 



